Solution to the Landau-Zener problem via Susskind-Glogower operators 
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We show that, by means of a right-unitary transformation, the fully quantized Landau-Zener 
Hamiltonian in the weak-coupling regime may be solved by using known solutions from the standard 
Landau-Zener problem. In the strong-coupling regime, where the rotating wave approximation is 
not valid, we show that the quantized Landau-Zener Hamiltonian may be diagonalized in the atomic 
basis by means of a unitary transformation; hence allowing numerical solutions for the few photons 
regime via truncation. 
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INTRODUCTION 



The Landau-Zener (LZ) problem [1H3 consists of a 
two-level system whose parameters are varied so that an 
anti-crossing of energy levels occurs. The transition be- 
tween two energy states at an avoided level crossing in 
two- and multi-level systems is one of the few exactly 
solvable problems of time-dependent quantum evolution 
043 ■ Dynamics in atomic, molecular and mesoscopic 
systems can be described by the LZ process; see, for ex- 
ample, references within 0, @- Recently, experimental 
realizations of many-body generalizations of LZ dynam- 
ics have been shown with ultracold atoms [|| and theo- 
retically analyzed considering strongly correlated bosons 
under fast sweeps @. The use of an interacting BEC 
driven in a bichromatic optical lattice has also been pro- 
posed as a realization of many-body nonlinear LZ dy- 
namics [lC 
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In particular, there exists an approximate solution for 
a non-interacting many-body generalization of the LZ 
problem, including coupling to a quantized field, suggest- 
ing that many-body LZ physics can be profoundly dif- 
ferent from the single two-level system interacting with 
a classical field case [l2|. Motivated by these results, 
we analyze the LZ problem of the one two-level system 
coupled to a quantized field. First, we introduce a quan- 
tized LZ model and propose a realization in circuit quan- 
tum electrodynamics (circuit-QED) where weak- and 
strong-coupling regimes can be obtained; we also dis- 
cuss plausible realizations of a formal many-body gener- 
alization of the model. Next, we present a right- unitary- 
transformation scheme to diagonalize the proposed quan- 
tized LZ model in the field basis under weak-coupling and 
show that the exact time evolution of the system is di- 
rectly related to solutions of the standard LZ problem. 
Finally, we show a parity-rclated-transformation scheme 
to diagonalize the proposed quantized LZ model in the 
two-level system basis under strong-coupling and show 
that the resulting infinite set of differential equations is 
amenable for numerical solutions for a few starting exci- 
tations in the field. 



MATHEMATICAL MODEL AND PHYSICAL 
REALIZATION 

A two-level system, that is, a qubit, driven by a quan- 
tized field is described by the model Hamiltonian, 

Hdf = -^qVz + FiLoa) a + hg(a + a^)cf x . (1) 

Where the Pauli matrices associated to the qubit are 
given by the operators &j, with j = z,±, and the sym- 
bol a (cv) is the creation (annihilation) operator of the 
quantized field. The qubit two-level transition and the 
field frequencies are given by loq and u, respectively, and 
the qubit-field coupling by the parameter g. 

In the weak-coupling regime, where the values of the 
qubit-field coupling are at least an order of magnitude 
less than the qubit transition frequency, g < O.Iloq, the 
rotating wave approximation (RWA) is valid and the 
Jaynes-Cummings (JC) model [l^ . in units of h, de- 
scribes the system, 



fr w „ - 1 - i \ 
ti.ic = ~^~ cr z + wa 1 o + A 



a' a- 



(2) 



The excitation number ./V = a T a + a z /2 is conserved by 
the JC model, [N, Hjc] = 0. Hence in a frame defined by 
the conserved excitation number rotating at the field fre- 
quency, to, by considering a time independent coupling, 
A, and a time dependent detuning between the qubit and 
the field frequency given by the original Landau-Zener 
(LZ) process, A = luq — to = —tluq, it is possible to write 
a quantized LZ, or LZ-JC Hamiltonian, in units of Hljo, 



H = —t <7, 



g (a'(T_ 



+ aa + ) 



(3) 



Where a scaled time r = v 2 t and according coupling 
g = X/v 2 has been used. This particular choice of pa- 
rameters means that the driving field is detuned to the 
blue of the qubit transition, w = (1 + t)oj = (1 +v 2 t)uJo; 
the positive constant v 2 is just the steepness of the linear 
frequency ramp. Choosing red detuning instead of blue 
is equivalent to replace r — > —t in Eq.Q. The cross- 
ing of the spectra is given at r — > where resonance is 
reached as ljq = tof. Notice that the quantized field fre- 
quency may be kept constant while adequately detuning 
the qubit transition. 
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A superconducting qubit coupled to a strip line res- 
onator [l4-16| may be described by the quantized LZ 
Hamiltonian in Eq.© when the transition frequency of 
the charge (flux) qubit is varied linearly in time by a 
driving charge (magnetic flux) leading to loq = (1 — t)lo. 
Furthermore, a circuit-QED implementation may allow 
the strong coupling needed to go beyond the RWA fl~7j ]. 
For simplicity, our analysis of the strong coupling case 
will use the Hamiltonian, in units of hu, 



H c = f a z + a) a + g(a + a))d x 



(4) 



where the transition frequency of the qubit is tuned by an 
external charge (magnetic field) to vary linearly in time 
as f = u 2 t; again, the positive constant u 2 is just the 
steepness of the driving ramp. In this model, the blue 
(red) detunning between the field and the qubit is pro- 
vided by r < 1 (r > 1), and resonant driving is obtained 
at t = 1. 

A BEC in an asymmetric double well trap fli| . where 
the depth of one of the wells is made to vary linearly 
in time, may be described by a nonlinear version of the 
standard two-level Landau-Zener process. In the mean- 
field approximation, classical dynamics of nonlinear LZ 
tunneling in Bose-Einstein condensates (BEC) has been 
studied and separated tunne ling defined by the fixed 
points of the system found [19(. A two- modes BEC 
driven by a quantized microwave field delivers a quantum 
version of the condensate in an asymmetric double-well 
model [2(| 21 1; the modes may be defined by two hy- 



perfinc structures of the ground state of a given atomic 
species, the time dependent detuning is given by the de- 
tuning between the qubit transition frequency and the 
driving field, and the nonlinearity may be tuned down 
by Feschbach resonances and neglected. Such consider- 
ations deliver what we will call a Landau-Zener-Dickc 
(LZD) Hamiltonian, 



H 



LZD 



-t S z + gN~ 1 / 2 [cJS-+aS 



Where the angular momentum basis related to angular 
momentum operators Sj, with j = z, ±, describes the en- 
semble containing N q qubits. Also, the model in Eq.© 
may describe the effective motion of a laser driven con- 
densate, under the two-mode approximation, coupled to 
an optical cavity 22H24j. when the detuning between 
the laser pump and cavity frequencies with respect to 
the difference between the energies of the two-lowest- 
momentum modes is set to vary linearly in time. While 
an approximate solution for the LZ problem is known for 
the LZD model 12} . an exact solution for its time evolu- 
tion is feasible but, at least for the time being, we restrict 
our analysis to the models of the one qubit within and 
without the rotating wave approximation described by 
EqQ and EqQ, respectively. 



DIAGONALIZATION IN THE FIELD BASIS 

In order to provide a time evolution operator for the 
quantum Landau-Zener Hamiltonian within the rotating 
wave approximation, Eq.Q, we diagonalize it in the field 
basis by using the Susskind-Glogower [25[ operators, 



V 



a T a+ 1 



-1/2 



a, 



(0) 



which arc rig ht-unitary, that is, VV^ = 1 and V^V = 
1 — |0)(0|. Their action on a Fock state of the number 
basis is given by V\n) = \n — 1) and V^^n) = \n + 1). 

Via the Susskind-Glogower operators, the LZ-JC 
Hamiltonian in Eq. (J3J) may be re- written as, 



H = tH + f^H LZ f, 



(7) 



where the auxiliary, Hq, and LZ-like, Hlz, Hamiltonians 
are 



F = (I-^)|0)(0|/2 

1 /2 

Hlz = -t a z + g (a^a + l) a x 



(8) 
(9) 



and the right-unitary transformation T is defined by [2 



f= I - a z + (I - <7 Z ) V /2 



(10) 



Notice that no approximation has been made in re- 
writing Eq.Q as Eq. ([?])• If a semi-classical quantization 
of the field were followed, like that proposed in Ref.(29j 
for time independent JC and Dickc models and only valid 
for coherent staes of the field, an approximate model lack- 
ing the Hq term would be obtained. 

As the auxiliary Hamiltonian Hq commutes with 
the transformed LZ-like Hamiltonian, T^H^zT, at 



any given time, 



nHo^HLz^f 



0, and 



using the fact that f tfl" LZ ( T i)f, ftff LZ (T 2 )T 



(5) ft 



T, it is possible to write the time 



evolution operator of the system described by Eq.© as 
U(t) = U (T)f^U LZ (T)f; (11) 
where it is trivial to find U {t) = e~ iT f ^° dt = e ~ iT ' 2Ha / 2 



due to the fact that 



t 1 Hq,t 2 H c 



= 0. 



Notice that the right-unitary transformation T acting 
on the dressed state basis yields, 

f(c n .o\n+ 1,0) ±c n .i\n, 1)) = \n) (c„, |0) ±c„,i|l)), 

(12) 

where the shorthand notation \n,x) = \ri)\x) = \n) fi e id® 
\x)atom with n = 0,1,2,... and x = 0,1, has been 
used. The time evolution operator for the LZ-like pro- 
cess, Ulz(t), in the dressed state basis is given by the 
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following set of coupled differential equations, 

1/2 

i d T c nt i{T) + tc, iA (t) -g(n+l) c„ j0 (r) = 0, 

1 /2 

i d T c nfi {T) - TCnfiir) - g(n + 1) c Mj i(t) = 0, 

(13) 

where the shorthand notation d T - denotes the partial 
derivative with respect to the scaled time r. Notice that 
this differential system is equivalent to that given by the 
standard two- level Landau-Zener process [4|, with the 
difference that the coupling between the two levels, g, 

1/2 

is enhanced by a factor (n + 1) due to the quantized 
field. 

The system of coupled differential equations separates 
into, 



[d 2 T + r 2 +g 2 (n + l)-i] c„ a (r) = 0, 
[d 2 + t 2 + g 2 (n + 1) + i] c„, (r) = 0, 



(14) 



which might be reduced to well-known differential equa- 
tions accepting Whittakcr functions, parabolic cylinder 
functions, and confluent hypergeometric functions of the 
first kind as solutions [3CH33| . The latter will be preferred 
for the sake of simplicity at t = 0; that is, even solutions 
evaluate to a non-zero constant while odd evaluate to 
zero. The properties of the hypergeometric function of 
the first kind, ii ? i(-, •, •), yield the time evolution for the 
amplitudes, up to a normalization and initial conditions 
factor, 



Cl , e (r) = e- <r 2/2 Q + 

c i,o(r) 



1 , i<£ 1 ■ 
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c , e (r)=e- 2 / 2 ^(f \,ir 2 ) 



C0,o( T ) 



-iT 2 /2 TP I 1 , 3 . ; 

1 1 ' 2 ~4~'2' 



(15) 



where the functions imply time and photon number de- 
pendence, Ca jP (r) = c x ^ p {g,n,T) with x = 0, 1 and 
p = e,o, the auxiliary characteristic values are defined 
as g n = g(n + l) 1 ' 2 . The time evolution operator Ulz is 
defined by the matrix elements, 



U 



(16) 



with 



Ul,l 
"1,2 
"2,1 
«2.2 



CO,e( T o) c l,e( T ) ~ c 0,o(t )ci ;O (t), 

ci, e (To)ci j0 (r) - ci )0 (ro)ci,e(r), 

Co,e(To)c ,o(T) - Co i0 (To)co,e(r), 

ci, e (ro)co,e(r) - ci (ro)co )0 (r), 



(17) 
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Figure 1. (Color online) Exact evolution of the population 
difference (<Tz), Eq.{T5]), for interactions starting at time (a) 
to — > — oo and (b) to = — 10 under weak coupling g = O.li^o 
and initial state \ipg) = \n,g > with n = 1 (solid black), 
n = 11 (dashed blue), n = 31 (dotted red) and n — 101 
(dot-dashed green). 



where 



is given by 



(t,to,gn) and the normalization factor 



7 = [cqA t o)ciA t o) ~ co )0 (to)ci )0 (t )] 



(18) 



with the number operator defined by n = era. The time 
evolution operator Ulz is exact, no approximation has 
been done. 

As mentioned before, the presence of the quantized 
field enhances the qubit-field interaction simulating a 
standard two-level LZ process with effective coupling g^; 
this can be observed graphically from the time evolution 
of the population difference, 



(a z {r)) = 0j(t )\U\t)<t z U(t)4;(t q )), 

= ^(r )\f^Ul z (r)a z UL Z (T)f\^(r )) 



(19) 



For the sake of historical comparison, FigUJa) shows 
the time evolution of the population difference, (ct z (t)), 
in the context of the original LZ process where the in- 
teraction starts at Tj — > — oo. Initial states given by 
\ip{n -> -oo)) = \n,g) withne {|1), |11), |31), |101)} are 
considered. The results are exact numerics from Eq.(fT9|) 
by using Eq. fTBT) . Figure [IJb) shows the finite time effect 
described by the exact time evolution. Initial conditions 
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Figure 2. (Color online) Asymptotic probability to find the 
qubit in the excited state, P e (j —¥ oo), for the symmet- 
ric crossing defined by a starting (ending) time to — > — oo 
(t — > oo) and the asymmetric crossing defined by tq — > — fO 
(r — >• oo) under weak coupling g — 0.1o;o and initial state 
jf/>o) = \n,g >. Plots are calculated both with the analytical 
asymptotic expansion given by Eq. (|20[> and from numerics 
of the exact evolution, Eq. (|19p by substituting the hypergeo- 
metric function by its asymptotic expansion up to third order, 



with initial times to — ► 
r -> 10 6 (Plots overlap). 



-10° and t 



TO and final time 



are the same described above, but for an initial time of 
interaction To = — 10. As expected jij], the population 
difference oscillates as soon as the finite time interaction 
starts. 

Figure [2] shows the asymptotic behavior for the ex- 
cited state probability, P e , for the symmetric crossing 
with starting (ending) times tq —> — oo (r — > oo) and the 
asymmetric crossing given by tq — > — 10 (r — > oo). The 
qubit is initially taken in the ground state and the field 
at an arbitrary Fock state. The probabilities are calcu- 
lated from the asymptotic expansion of Eq. (fT!?|) . via the 
parity and asymptotic properties of the hypergeometric 
function 31|, and both yield the expression 



P e (r oo) = (1 + (<j z (t -> oo)))/2, 



1 



(20) 



which is equivalent to the tunneling probability predicted 
in the standard LZ process. The tunneling probability is 
enhanced by the number of photons in the quantized ficd 
as g n =g(n+l) 1 / 2 . 

The exact time evolution operator found in this section 
describes any given set of parameters involving initial 
system state, weak coupling, initial and final time. Fig- 
ure [3] shows the time evolution of the population differ- 
ence, (o- z (t)), in a purely asymmetric case starting from 
the crossing. The qubit is initialized in the excited state 
|e) and the field in the Fock states n G {0, 10,30, 100}]; 
short and long interaction times are shown in subfigurcs 
FigEI a ) and FigJ3Jb), respectively. In the asymptotic 
infinite time, the probability of finding the qubit in the 
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Figure 3. (Color online) Exact evolution of the population 
difference (<j z ), Eq. (|19|I . for interactions starting at time to = 
under weak coupling g — O.lt^o and initial state |^o) = 
|n, e > with n = (solid black), n = 10 (dashed blue), n = 30 
(dotted red) and n — 100 (dot-dashed green) for (a) short and 
(b) long interaction times 



1 ± e^ 9 "' 2 ) /2, (21) 



excited state is 

Pe{n =0,Tf -> OO) 

depending if the qubit starts in the excited (plus sign) 
or the ground state (minus sing). This probability is is 
plotted in Fig0]for the qubit starting in both the excited 
and the ground state. 

DIAGONALIZATION IN THE ATOMIC BASIS 

For moderate and strong coupling, g > O.lwo, the ro- 
tating wave approximation is not valid and the complete 
Hamiltonian in Eq.([3]) has to be considered. Via the uni- 
tary transformation, 



i? 



D — i-n-(iT x — 1)0^0/2 



(22) 



the complete Hamiltonian in Eq.Q becomes diagonal in 
the qubit basis, 

R^H C R 



— t z cos(7ra'a) + a} a + g(a + a T ); (23) 



H, 



The use of the unitary transformation in Eq. (|22[) is equiv- 
alent to consider a parity chain basis and the corre- 
sponding creation/annihilation operators as proposed in 
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Figure 4. (Color online) Asymptotic probability to find the 
qubit in the excited state, P e (r —¥ oo), for the interaction de- 
fined in FigOstarting from initial states given by |-0o) = \n, e) 
(solid line) and |^o) = \n, g) (dashed line). Plots are calcu- 
lated both with the analytical asymptotic expansion given by 
Eq, (|21[l and from numerics of the exact evolution, Eq. (|19[) by 
substituting the hypergeometric function by its asymptotic 
expansion up to third order, with r — > 10 6 (Plots overlap). 

Ref.[34| for the time independent qubit-field interaction 
in the strong-coupling regime. It is straightforward to 
see the relation between these two approaches from 

= aa x , (24) 

= tfa x , (25) 
Wa z R = & z cos(7ra^a) + ia y sin(7ra' a) 
= a z cos(7ra^a) 

= (-1)***., (26) 

as the parity operator, II = (— l) ata tr z , in the basis de- 
fined by the unitary Eq. (|22p is given by 

R^UR = (-I) a1 ' a i?V z i?, 

= (r z . (27) 

That is, the complete Hamiltonian Eq. ([2lt)) conserves par- 
ity, [ttfllR, R)H C R] = [a z ,H R ] = 0. 

By moving into the rotating frame defined by the free 
field, Uf{t) = e~ m aT , the dynamics are given by the 
Hamiltonian, Hrf = Hrq + Hm, with 

Hro = t a z cos(7ra^a) 

H m = g(ae lT + aV iT ). (28) 

The first of these terms, Huq, is diagonal in both the 
qubit and Fock basis and commutes with itself at dif- 
ferent scaled times, {Hr (ti), Hr (t2)] = 0; that is, it is 
possible to use a unitary transformation, 

U (t) = e-' lT ^~~ & V 2 : (29) 



such that the system is described by, 

tiRFU = g(p,e +a'e )e - y '. (30) 

This time dependent Hamiltonian produces two infinite 
sets of coupled first order differential equations for the 
field, one for each qubit state x G {0, 1}, 

i d T c xfi (T)=ge- lT e^ 2 CxA ( T ), (31) 

for n = and 

i d T c x , n (T) = 5 n 1 / 2 e ir e ±iT2 (- 1 )"" 1 c x , n _ 1 (r) 

+ g (n+l)^e— e ±^(-^ 1 Cxn+l{T) 

(32) 

for n > 1; for the sake of simplicity, dimension has 
been set to units fhui. The notation \<f>) = R\ip) = 
J2 X n c x,u(t)\x, n) has been used. For Fock states with 
photon number m, the differential set defined by Eq.([32|) 
may be truncated at an arbitrary large n ^> m, 

i d T c Xifi (r) = n^e^-^c^ir). (33) 

This is particularly helpful for initial states with small 
number of photons, in these cases numerical solutions 
may be given. Figure [5] shows numerics for the popula- 
tion difference, 

<a 2 (r)> = (ip(T)\a z cos (irtfa) |^(t)>. (34) 

An initial state \ipo(r = I)) = |0, e) is taken and the set 
of coupled differential equations is truncated at length 
one hundred, that is, m = 100. Qubit-field couplings in 
the range g e {0.1, 1, 3, 10}w are considered. The initial 
time t = I is chosen to emulate the case pictured in 
FigEl both systems start from resonance, loq = oj. The 
black solid line in FigJ3] and FigfS] represents identical 
initial conditions, |^o) = 0, e) and coupling g = O.Ilu 
and deliver similar dynamics. The rest of the couplings 
treated in FigJS] show dynamics similar to those under 
the RWA, Eq. ([7j), for small normalized times, r <C I 
and, then, the coherent oscillations break due the action 
of the counter-rotating terms, as expected. 

CONCLUSION 

We have presented a right-unitary approach to solve 
the qubit-quantized-field interaction under the rotating 
wave approximation with frequency detuning varying 
linearly in time, Eq.Q. The model may be realized 
in circuit-QED We have diagonalized the model 

Hamiltonian, Eq.Q, in the quantized field basis and 
shown that the procedure to obtain the non-trivial in- 
gredient of the evolution, U~Lz{t) in Eq. ((TT| . is already 
known from the interaction of a classical field with a 
qubit Q. The presented solution is exact. Its analytical 
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Figure 5. (Color online) Numerical time evolution of the pop- 
ulation difference (<t z ), Eq.lJMJ), for interactions starting at 
time to = 1 under couplings g — O.Iuj (solid black), g = lw 
(dashed blue), $ = 3a; (dotted red) and g = Wuj (dot-dashed 
green) and initial state |^o) = |0, e >. For the sake of compar- 
ison, the solid black line in this figure corresponds to initial 
conditions and parameters identical to the solid black line in 
FigH 



closed form allows its use in modular scenarios to engi- 
neer particular states or Hamiltonians; compare, for ex- 
ample, with Ref.[3|, where different symmetric and asym- 
metric crossings in the standard Landau-Zcner model 
are pro posed and may be used for state engineering, or 
Ref. [351] . where modular cavity-QED is proposed to engi- 
neer exotic lattice systems. In the asymptotic symmetric 
case equivalent to the standard LZ problem, the quan- 
tized version for initial separable states presents a similar 
asymptotic behavior for the probability of the LZ tran- 
sition, P e (r -> oo) « 1 - e^ 9 ' in Eq.flSOJl, with the 
distinction of an enhanced coupling proportional to the 
square root of the number of photons in the initial state, 
g A = g(n + l) 1 / 2 . 

The strong coupling dynamics of the system, where 
the rotating wave approximation is not valid, has been 
studied via an unitary transformation that diagonalizcs 
the Hamiltonian Eq.((4]) in the qubit basis. This operator 
approach is comparable to defining a parity chain basis 
for the system [34|. The system dynamics is given by 
an infinite set of coupled differential equations amenable 
to numerical solutions, via truncation, for small initial 
number of excitations in the quantized field. 
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